Bose-Einstein condensates with attractive interactions on a ring 
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Considering an effectively attractive quasi-one-dimensional Bose-Einstein condensate of atoms con- 
fined in a toroidal trap, we find that the system undergoes a phase transition from a uniform to a 
localized state, as the magnitude of the coupling constant increases. Both the mean-held approxi- 
mation, as well as a diagonalization scheme are used to attack the problem. 
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Effectively attractive Bose-Einstein condensates have 
concentrated a lot of attention recently, in connection 
with the experimental formation of bright solitons in 
them [1,2]. More specifically, in the experiments of Refs. 
[1,2], 7 Li atoms were confined in quasi-one-dimensional 
traps. With use of the Feschbach resonances [3] the ef- 
fective coupling constant that describes the atom inter- 
actions was then tuned and as it became negative - cor- 
responding to an effective attraction between the atoms - 
localized states, "bright soliton trains" were observed to 
form in Ref. [1] , while a single bright soliton was observed 
in Ref. [2]. References [4,5] have examined theoretically 
these systems. In the limit where transversely to the 
long axis of the trap the gas is in the lowest harmonic- 
oscillator level, the transverse degrees of freedom are 
frozen out and the system is essentially one-dimensional 

Motivated by these developments, we study here an 
effectively attractive one-dimensional Bose-Einstein con- 
densate of atoms confined in a toroidal trap. In a recent 
theoretical paper Kanamoto, Saito, and Ueda have inves- 
tigated the ground state and the low-lying excited states 
of such a system [5] (see also Ref. [7] for a detailed dis- 
cussion of this problem). As shown, the Gross-Pitaevskii 
mean-field theory predicts a quantum phase transition 
between a uniform state and a localized state as the ab- 
solute value of the strength of the interaction inreases. 
Furthermore, numerical diagonalization of the Hamilto- 
nian for a finite number of bosons shows that the tran- 
sition in this case is smeared out, as expected in finite 
systems. 

In our study we examine the same problem using dif- 
ferent techniques. Initially we use the mean-field approx- 
imation with a properly chosen variational wavefunction 
to study the phase transition and the order parameter 
in the two phases. Then, working in the same truncated 
space we use a Bogoliubov transformation to diagonalizc 
the Hamiltonian. Having diagonalized the problem, we 
examine the lowest state of the system, the low-lying ex- 
cited states, as well as the depletion of the condensate at 
the region of the transition. Our results are consistent 
with those of Ref. [5] . 

Let us therefore consider a Bose-Einstein condensate in 
a toroidal trap. Following Ref. [5], we assume that the 
system contains A bosons, that the radius of the torus is 



R and its cross section is S = irr 2 , with r <C R. If i>(9) 
is the field operator, the Hamiltonian is 



H = 



ffl 



(i) 



where Uq = 8ttciR/S, with a being the scattering length 
for elastic atom-atom collisions, and 9 is the azimuthal 
angle. Here the length is measured in units of R and 
the energy in units of Ti 2 /2mR 2 , with m being the atom 
mass. 

Let us start with the mean-field approach. Within this 
approximation the system is described by a single wave- 
function, the order parameter i/j(t), and the many-body 
state is the product nj^(rj), with i = 1, . . . , A, where 
N is the number of atoms. Therefore this approximation 
ignores correlations between the atoms and in general it 
has a higher energy than the exact solution that one can 
get by diagonalizing the Hamiltonian. 

In this problem it is natural to work in the basis of 
plane-wave states <pi(9) = e tW /V2tt and according to 
the analysis of Ref. [5] , the order parameter close to the 
transition consists of the state (fro (the dominant compo- 
nent), and the states <j>±\. To get a simple physical pic- 
ture, we thus develop a variational approach, expanding 
ip(9) in the basis of the (fri states and keeping only these 
three components. This is a reasonable assumption, since 
states with higher values of I have higher kinetic energy, 
and indeed as shown in Ref. [5] this approximation gives 
qualitatively (but not quantitatively) good results. We 
thus write 



%p(9) = c_i0_i + c (fr + c\§\ 



(2) 



Since the dominant component of ip is <f>o, therefore 
| c 1 ^> |c_i|, and |c | ^> \ci\. Because of the symmetry of 
the problem, |c_i| = |ci|, which also guarantees that the 
total angular momentum of tfr is zero, as it should. The 
normalization condition imposes the further constraint 
Ic^xl 2 + |c | 2 + | Cl | 2 = 1. 

To proceed, we express the energy per particle to m 
terms of the coefficients Ci, which are then determined 
by minimizing £q with respect to them [8], 
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c = 2| Cl | 2 + |[|c | 4 + | C _!| 4 + \c\\ + 4|c | 2 | C _ 1 | 2 
+4|c | 2 | Cl | 2 + 4|c_ 1 | 2 |c 1 | 2 + 2clc*_ lC l + 2(c 2 )*c_ lCl ], (3) 

where 7 = NU /(2ir) — ANaR/S is essentially the ratio 
between the interaction energy and the kinetic energy. If 
Cj = \cj\e : ' , then the last two terms in Eq. (3) are equal 
to 4|co| 2 |ci| 2 cos(#_i + 81 — 26* ). To minimize e one has 
to choose 6L1 + 6»i - 29 = 0. Since |c | 2 = 1 - 2|ci| 2 , 



; (l + 2 7 )-7 7 | Cl | 



(4) 



Here 7/2 is the interaction energy per particle of the 
uniform state. The above equation implies that for 7 > 
7 cr , where 7 cr = —1/2, the minimum of the energy occurs 
for c = 1, and c_i = C\ = 0. In this regime the density 
of the system is uniform and the energy per particle is 
7/2. On the other hand, if 7 < 7 cr , e (|ci| 2 ) is like a 
mexican hat, and its minimum occurs for all three Cj 7^ 0, 
the energy per particle is lower than 7/2 and the cloud 
develops a non-uniform density. Our approach implies 
that the transition takes place for the same critical value 
of 7 as in the exact solution of the mean-field approach 
[5,7]. It is interesting that within our scheme including 
more basis states in the order parameter does not affect 
7 cr . More precisely, including the states with |Z| < m, 
the quadratic terms in the energy have the form 



eo - \ = 2| Cl | 2 (l + 2 7 ) + 2|c 2 | 2 (4 + 2 7 ) 
+ ... + 2|c m | 2 (m 2 + 2 7 ). 
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FIG. 1. The magnitude of the trial order parameter tp, 
Eq. (2) , calculated within the mean- field approach, with the 
coefficients \a\ given by Eq. (6). The phases of a are chosen 
to be equal. One phase is arbitrary, the other one is chosen so 
that the energy is minimum, and the third one is free (reflect- 
ing the rotational invariance of the problem) and it has been 
chosen so that the maximum of the wavefunction is at 6 = 0. 
The straight line corresponds to a uniform condensate, with 
7 > 7cr, the less localized state to 7 = —0.51, and the more 
localized state to 7 = —0.55. 



and we will calculate terms in the energy of order 1/N, 
let us expand the field operator -0, 



(5) 



ij){9) = <p_x c_i + 4>o c + 0i £1, 



(8) 



From the above expression it is obvious that 7 cr is always 
— 1/2, essentially because of the high kinetic energy of 
atoms in states with high values of I, that scales as I 2 . 

Let us now examine the non-uniform state when 7 < 
7 cr . Demanding that deo/d\ci\ 2 = 0, we get 

2 4 
|c±i| 2 = 7~(T _ |c | 2 = 1 - - 7cr)- (6) 

Figure 1 shows the corresponding order parameter. In 
addition 



7 4 



7cr) 2 . 



(7) 



If 7 = —1/2 — S, where S — > and S > 0, the above 
equation implies that e = —1/4 — 5/2 — 8<5 2 /7, which 
should be compared with the exact mean-field solution 
[5], Cq x ~ —1/4 — 5/2 — 2S 2 . Therefore our trial wave- 
function reproduces the linear correction, but it fails in 
the quadratic term and it has higher energy. This is not 
a surprise, however, since our approach is variational and 
it does not necessarily capture all these features. Includ- 
ing more basis states, one can improve the wavefunction 
and the energy to the desired order in S. 

Turning to our alternative approach of the problem, 
where we will go beyond the mean-field approximation 



where q is now the annihilation operator of an atom with 
angular momentum I. Again we have restricted ourselves 
to the states with I = 0, ±1. In this basis the Hamiltonian 
(1) can be written as 

H = ^ + 4^ ^2 Ck4crnC n 5 m+n -k-l- (9) 



klmn 

A similar Hamiltonian has been studied in Ref. [8] in the 
context of weakly-interacting Bose-Einstein condensates 
under rotation. Denoting the basis vectors as 

\m) = IN-uNo,^) = |(-l) m ,0 Ar - 2m ,(+l) m ), (10) 

where TV; is the occupancy of the state <f>i, with J^i Ni = 
N, and ZiVj = 0, it is straightforward to calculate the 
diagonal matrix elements, 

(m\H\m) = 2m+—[(N- 2m)(N -2m - 1) 

47T 

+2m(m- 1) + 8m(N - 2m) + 4m 2 }. (11) 

Assuming that the system is close to the transition point, 
but on the side of the "uniform state" (in the infinite- A'' 
limit of the mean- field approximation), then m is of order 
unity and thus m <C N. In this limit, 
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(m\H\m) w ±{N - 1) + 2m(l + 7), 



(12) 



where j(N — l)/2 is the total energy of a condensate 
with uniform density. In Eq. (12) terms of order m 2 have 
been neglected, and therefore this approach is not valid 
for 7 < 7 cr (since m becomes of order N there). 
Turning to the off-diagonal matrix elements, 

Uo. 



(m\H\m + 1) = -^2y/(N - 2m){N - 2m - l)(m + l) 2 , 

(13) 

and in the limit m -C N, 

(m\H\m+l) =7(m + l). (14) 
From Eqs. (12) and (14) can be written as 

£-l(iV-l) = 

(1 +7)(c t _ 1 c_i +cjci) + t(c-ic! +c_ici). (15) 

The above Hamiltonian can be diagonalized with use of 
a Bogoliubov transformation [9], 



b = Aic_ a + A2C1 
d = A 2 c_i + Aic{. 



(16) 



For and d to satisfy bosonic commutation relations, 
A 2 — A 2 = 1. Following the usual tricks, i? can be written 
in the diagonal form 



H - ^(JV -1) = E + E(tfb + $d), (17) 



where 



and 



S (7) = 2[(l+7)A?-7AiA 2 ], (18) 



E{ 1 ) = {l+ 1 )(X 2 1 +X 2 2 )-2 1 X 1 X 2 . (19) 
In addition, to eliminate the off-diagonal terms one has 
2AiA 2 7 



A? + A 2 7 + 1 ' 



(20) 



Parametrizing, i.e, writing Xi = sinh(6>) and A 2 
cosh(#), Eq. (20) can be written as 



tanh(20) =7/(7+1), 



(21) 



and in order for a solution to exist, 7 > —1/2. Equation 
(21) implies that 9 = ln(27 + 1), which can be used to 
express all the quantities in terms of 7. Thus we get for 
the eigenvalues of H, 



£n b ,n d ~ j(N - 1) = -(7 + 1) + Vl + 2 7 (1 +n b + n d ), 

(22) 



where nt,, nd = 0, 1, 2, . . . are the eigenvalues of the num- 
ber operators b'b and d)d respectively. For N — > 00, 
£0,0 = 7/2 = e , in agreement with mean field. 

Having diagonalized H, we can easily get the occu- 
pancy of the 4>i=±i states, which is 



(4iC±i)/^V = \\/N. 
Expanding we get that close to the transition, 

(4ic±i) ^ V2( 7 + l/2)- 1 / 2 



N 



N 



(23) 



(24) 



Since the occupancy of the <f>±\ states is a smooth and 
continuous function for any value of 7 when N is finite, 
and since in the limit N — > 00, (c± 1 c±\)/N has to equal 
l c ±i| 2 [10], a simple geometric construction shows that 
the value of the function in Eq. (24) at 7+ = —1/2 + 5, 
minus the value of |c±i| 2 at 7_ = —1/2 — 5 has to be 
proportional to 8, and thus 



(25) 



in agreement with Ref. [5], where it was found numeri- 
cally that n w 1.077. Therefore to leading order 



and 



|c±i| 2 = - 7 (7 + l/2 + ^- 2/3 ), 



7cr = -1/2 - V N- 2 / 3 . 



(26) 



(27) 



Turning to the energy spectrum, in the lowest state 
fib = rid = 0, and thus the ground state energy is 



£o(7) = -(7 + 1) + (27 + 1 



,1/2 



(28) 



measured from the energy of the state with uniform den- 
sity, j(N-l) /2. Making use of Eq. (27) we plot the func- 
tion £0(7 + 7]N- 2 / 3 ) in Fig. 2 for TV = 100 (top curve), 
500 (middle curve), and N — > 00 (bottom curve), with 
T] = 1.077. From Eq. (27) we get that to leading order 
£ ( 7 = -1/2) = -1/2 + (2r 7 ) 1 / 2 7V-i/3. 

The low-lying excited states of H lie above the lowest 
state by E(-f)(rib + rid), where 



£( 7 ) = (2 7 +l) 1/2 . 



(29) 



The energy of the lowest excited state is thus E(-f), in 
agreement with Ref. [5]. Figure 3 shows E(j + -qN~ 2 / 3 ) 
for TV = 100 (top curve), 500 (middle curve), and 
N -> 00 (bottom curve). From Eq. (27), £(7 = -1/2) = 
(2r ] ) 1 / 2 N- 1 / 3 . 

Finally the depletion of the condensate AN is, 



AiV = (c^c-i) + (clci) = 2\f - ( 7+ 1 



VV27TT 



1 • (30) 



Figure 4 shows A7V(7 + t]N^ 2/3 ) for N = 100 (bottom 
curve), 500 (middle curve), and TV — > 00 (top curve). 
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According to Eq. (27), to leading order AA(7 = -1/2) = 

(x/274)7T 1/2 AT 1/3 - 

To summarize, we examined a Bose-Einstein conden- 
sate in a toroidal trap, with an effective attractive inter- 
action between the atoms. Using a variational approach 
and working within an appropriately chosen set of basis 
states, we demonstrated that within the mean-field ap- 
proximation for a strong enough coupling constant a con- 
densate of uniform density becomes unstable against the 
formation of a localized state. Going beyond the mean- 
field approximation, we diagonalized the Hamiltonian in 
the same truncated space just above the transition and 
we calculated the energy of the low-lying states, as well as 
the depletion of the condensate. While the energy in this 
case is lower than that of the mean-field, the two meth- 
ods give the same results in the N — > oo limit, assumed 
implicitly in mean-field. 
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FIG. 2. The energy E (-y + ■qN' 2/ ' i ) of the lowest state of 
H, Eq. (28), for N = 100 (top curve), 500 (middle curve), and 
N — > oo (bottom curve), with rj = 1.077. 
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FIG. 3. The energy £(7 + t/A~ 2/3 ) of the first excited state 
of H, Eq. (29), for N = 100 (top curve), 500 (middle curve), 
and N — ^ oo (bottom curve), with r\ — 1.077. 
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FIG. 4. The depletion AiV(7 + r]N' 2/a ) of the condensate, 
Eq. (30), for N = 100 (bottom curve), 500 (middle curve), 
and N — > oo (top curve), with r\ — 1.077. 
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